Coherent response of the Hodgkin-Huxley neuron in the high-input regime 
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We analyze the response of the Hodgkin-Huxley neuron to a large number of uncorrelated stochas- 
tic inhibitory and excitatory post-synaptic spike trains. In order to clarify the various mechanisms 
responsible for noise-induced spike triggering we examine the model in its silent regime. We report 
the coexistence of two distinct coherence resonances: the first one at low noise is due to the stimula- 
tion of correlated subthreshold oscillations; the second one at intermediate noise variances is instead 
related to the regularization of the emitted spike trains. 
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A. Introduction 

In the last decades a large number of studies have been 
devoted to the characterization of the response of simple 
and more elaborated neuronal models under the influ- 
ence of a large variety of stochastic inputs 0, l1 US ■ One 
of the main reasons justifying the interest of neurosci- 
cntists for this subject resides in the observation that 
in vivo neocortical neurons are subjected to a constant 
bombardment of excitatory and inhibitory post-synaptic 
potentials (EPSPs and IPSPs), somehow resembling a 
background noise Q. 

Among the many proposed biophysical models the 
Hodgkin-Huxley model [5J can still be considered as a 
reasonably valid framework for exploring the main fea- 
tures of neural dynamics |16|. In order to understand 
the origin of the variability observed in the distribution of 
spikes emitted by cortical neurons (2^ the response of the 
Hodgkin-Huxley (HH) model has recently been studied 
under the influence of additive noise [lj, l2Ct [25| and sub- 
jected to trains of post-synaptic potentials [JM H3 ■ 

One of the most interesting phenomena observed ex- 
perimentally and numerically for excitable neuronal sys- 
tems driven by noise is coherence resonance (CR) |cill7|: 
it refers to the regularization of noise- induced oscillations 
occuring for finite (not zero) noise intensity in the ab- 
sence of an y ex ternal injected signal (for a comprehensive 
review see pj|). Evidences of CR have been reported ex- 
perimentally for the cat's spinal and cortical neural en- 
sembles |l5j and theoretically for various neuronal mod- 
els ^l). In particular, CR has been observed for the HH 
model [lOl |27| , but these results mainly refer to additive 
continuous noise. 

Our aim is to perform a detailed analysis of the re- 
sponse of the HH model in the silent regime when it 
is subjected to many (hundreds or thousands) stochas- 
tic trains of post-synaptic potentials (PSPs) per emit- 
ted spike (i.e . the system is in the so-called high-input 
regime |3l l23jp . 

The different mechanisms responsible for noise- induced 



neuronal firing are analyzed in terms of statistical in- 
dicators (interspike-interval distributions and their first 
moments) as well as of dynamical indicators (autocorre- 
lation functions). The analysis of the correlation time 
reveals a CR associated to the stimulations of coherent 
subthreshold oscillations around the rest state. This new 
type of CR coexists with standard CR related to the reg- 
ularization of spike sequences [TH ]. 



B. Model and methods 

The HH model describes the dynamical evolution of 
the membrane potential V(t). It can be written as 

C^r = g Na m 3 h(E Na -V) + g K n 4 (E K -V) + (1) 



dt 



+ g L (E L -v) + i(t) 



where I{t) represents the synaptic current and the evo- 
lution of the gating variables X — to, n, h is ruled by 
three ordinary differential equations of the form dX/dt = 
ax(V)(l — X) — Px(V)X. The parameters entering in 
Eq. and the expressions of the nonlinear functions 
ax(V) and Px{V) are reported in [l2|. 

In this study we consider the response of the single HH 
model to N e (resp. Ni) uncorrelated trains of EPSPs 
(resp. IPSPs) with NJN) ranging from 10 to 10,000. 
By following Refs. 0,0, each PSP is schematized as an 
instantaneous variation of the membrane potential by a 
positive (resp. negative) amount AV for excitatory (resp. 
inhibitory) synapses. The amplitude of each voltage kick 
is assumed to be 0.5 mV, i.e. reasonably small (s=s 7 %) 
with respect to the distance between the "threshold" for 
spike initiation for rapid EPSPs and the resting potential 
(w 6-7 mV) estimated for the HH model [01 Hf]. More- 
over, amplitudes « 0.5 mV are comparable with average 
EPSPs experimentally measured for pyramidal neurons 
in the visual cortex of rats [|| . This amounts to exciting 
the neuron with an impulsive current 
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N e Ni 

m = q[EE^-4)-EE^- o] (2) 



fc=i 1 



m — 1 n 



2 



where t k (resp. t™J are the arrival times of the EPSPs 
(resp. IPSPs) and Q = CAV is the charge associated to 
each kick. In order to mimic the inputs received by corti- 
cal neurons, for each afferent (excitatory and inhibitory) 
synapse the time interval distribution between PSP in- 
puts is chosen to be Poissonian with an average frequency 
u = 100 Hz [H 1 . 

The stochastic input can be characterized in terms of 
the average and the variance of the net spike count within 
a temporal window AT 



A', 



7V(AT) = ^nf(AT)-^n^(AT) 



(3) 



k=l 



where n^(AT) (resp. n^(AT)) represents the number 
of afferent EPSPs (resp. IPSPs) received from neuron k 
(resp. m) in the interval AT. In the high input regime, 
by assuming uncorrelated input s pike trains, A^(AT) fol- 
lows a Gaussian distribution (cf. |12|). with average and 
variance derivable within the shot noise formalism as 



< JV(AT) >= (N e - N t )ATv a 



(4) 



and 



Var[N(AT)} = (N e + N^ATuq = a 2 ATv . (5) 



The parameter a = (N e + N) measures the standard 
deviation of the stochastic input process. 

The average current stimulating the neuron is given by 



/ = 



CAV < N(AT) > 
AT 



= CAVv (N e - Nt). 



(6) 



I represents the bifurcation parameter ruling the dy- 
namics of the deterministic HH model; for / < Isn — 
6.27/iA/cm 2 the model is in the silent regime, i.e., in the 
absence of noise the dynamics is attracted towards a sta- 
ble fixed point and the neuron does not fire. However, 
since the fixed point is a focus the relaxation towards the 
rest state occurs via damped oscillations (subthreshold 
oscillations) |19|. Periodic firing is observed for the HH 
model only above Isn and it is associated to frequencies 
in the range 50 — 80 Hz. 

We have verified that in the high input regime the re- 
sponse of the neuron (for fixed v$) is determined once 
fixed I and a, therefore it does not depend separately on 
N e and Ni, but only on their difference (N e — Ni) (x I 
and sum (N e + N) <x a 2 [T^ |. 

In order to characterize the output of the neuron and 
to examine the coherence effects in the response we have 
employed the following indicators: the distribution of the 
output interspike intervals Pisiit) and its first moments; 



1 It should be stressed that the model here analyzed is driven by 
impulsive post-synaptic currents and it should not be confused 
with conductance-driven models 20, 251. 



the coefficient of variation of the output interpike inter- 
vals distribution (Pisi(t)) defined as R = Sisi/Aisi, 
where Aisi is the average and Sisi the standard devia- 
tion of the Pisi(t) (for a perfectly periodic response we 
have R = and R = 1 for Poissonian output; the nor- 
malized autocorrelation function C(t) for the membrane 
potential and the correlation time [lj| defined as 



C 2 (t) alt 



CR is usually identified by a minimum (resp. a max- 
imum) in the coefficient of variation R (resp. in r c ) oc- 
curring at finite noise variance [Tl| . 



C. Results 

The HH neuron subjected to stochastic input exhibits 
a noise-induced (irregular) spiking behaviour with an av- 
erage firing rate v out = 1/Ajsi- For fixed I the firing 
rate increases with the standard deviation of the noise 
(cf. Fig.HJ). 

The noise-induced firing activity becomes more and 
more pronounced by approaching the repetetive firing bi- 
furcation Isn for the deterministic HH. The HH model 
is a type II neuron, therefore the onset of oscillation at 
/ ~ Isn is associated with a finite frequency (~ 50 Hz), 
however, in presence of noise, arbitrarily low spike rates 
can be observed even in the silent regime. 
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FIG. 1: Silent regime: average firing rate function 
of a for increasing values of I (from the bottom to the top 
/=-l,0,l,2,3,4,5,6.15 fiA/cm 2 ). 



1. Response of the HH model in the low and high noise 
limit 

Let us now examine more in detail the origin of noise- 
induced spiking in the low and high noise limits. For low 
noise we observe the coexistence of multi-peaks and an 
exponential tail in the Pisi{t) (cf. Fig. [SJ. As explained 
in the following the multi-peak structure is related to 
the relaxation dynamics of the membrane potential V(t) 



3 



following a spike emission, while the exponential tail is 
associated with noise induced activation processes. 



T 
/fp 



20 40 60 80 100 

t (ms) 



FIG. 2: Pisi(t) obtained for the stochastic input © with 
a = 4.6 (continuous curve) and the relaxation dynamics of 
the membrane potential V(t) following a step current stimu- 
lation (dashed line). The position of the output spike has been 
shifted to t = and the action potential has been rescaled to 
better reveal the relaxation oscillations. In the inset the ex- 
ponential tail due to the activation mechanism is shown. The 
data refer to / = 6.15/iA/cm 2 . 

To better understand the first mechanism we have 
studied the evolution of V(t) following a step current of 
amplitude I, i.e., 



I(t) = 



q if t < o 
/ if t > o 



(8) 



As shown in Fig[2] (dashed line) V{t) exhibits one spike 
followed by damped oscillations. This suggests that the 
probability of eliciting a second spike is enhanced in cor- 
respondence the maxima of the relaxation oscillations, 
where V(t) is nearest to the firing threshold. Moreover, 
it is possible to show that the period T n \ of the first os- 
cillation has a clear nonlinear origin, while the period of 
successives oscillations 7] can be explained by a linear 
stability analysis around the fixed point. Indeed we have 
verified that the first peak (7/ p ) observed in the distribu- 
tion Pisi(t) (continuous curve in Fig|2J) is related to T n i, 
while the successives peaks in the Prsrit) correspond to 
the linear oscillations of period 7] [l2|. 

The second firing mechanism responsible for the ex- 
ponential tail in the Pisi(t) is due to the competition 
between two effects, on one side the tendency of the dy- 
namics to relax towards its stable fixed point and on the 
other side noise fluctuations that instead lead the sys- 
tem towards the firing threshold. In this framework the 
dynamics of V(t) can be described as the overdamped 
dynamics of a particle in a potential well under the in- 
fluence of thermal fluctuations and the average firing 
time, or activation time Tg. can be expressed in terms of 
the Kramers expression |h4] for sufficiently small noise: 



T a tx e 



W/<j 2 



(9) 



where a 2 plays the role of an effective temperature and W 
of an energy barrier. For a 1 < W , i.e., in the low noise 



limit, the dynamics can be characterized as an activation 
process: Pisi(t) resembles a Poissonian distribution with 
R « 1. For a 2 > W, i.e., in the high noise limit, the 
multi-peak structure and the exponential tail disappear 
and Pisi(t) reduces to an inverse Gaussian corresponding 
to the distribution of the first passage times for a diffusive 
process plus drift |26|. 

As a further indicator we have estimated the spike trig- 
gered average potential (STAP) 0, that gives the aver- 
age shape of V{t) preceding the emission of a spike, for 
sufficiently long ISIs and for small a. As shown in Fig. [3 
V(t) oscillates with period 7] before firing; therefore 
the emission of a spike (for long ISIs) is triggered by the 
excitation of linear subthreshold oscillations around the 
fixed point. Thus the HH neuron nearby the rest state 
acts as a sort of selective filter since it responds (by emit- 
ting a spike) with higher probability when it is stimulated 
with a specific input frequency (w 1/2] = 61 — 88 Hz for 
< I < Isn)- This result agrees with the analysis of 
where it has been shown that a silent HH neuron sub- 
jected to a sinusoidal current optimally resonates when 
forced with a frequency linearly correlated with 1/T;. 
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FIG. 3: / = 5fiA/cm 2 and a = 5.7. STAP preceding the 
emission of a spike. At time t = the potential overcomes 
the spike detection threshold (fixed at —5 mV). STAP has 
been calculated by averaging only over ISIs longer than 150 
ms. 



2. 



Coherent response of the Hodgkin- Huxley neuron 



In the silent regime we have found the coexistence of 
two CRs: one is related to the regularization of the emit- 
ted spike trains at intermediate noise levels and has previ- 
ously been reported in fToll27|: the second one, observed 
for the first time, occurs at very low noise and is associ- 
ated to the coherence of subthreshold oscillations. 

The first CR can be explained by noticing that the 
dynamics of excitable systems can be characterized in 
terms of two characteristic times 0] : an activation time 
T a , representing the time needed to excite the system, 
and an excursion time T e , indicating the duration of the 
spike (i.e., the time spent in the excited state). An ISI is 
given by the sum of these two times and therefore R can 
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be (formally) expressed as the sum of two parts depend- 
ing separately on T a and T e . These two contributions 
vary in an opposite way when the noise is increased and 
the competition of these two tendencies leads to the max- 
imal coherence (associated to a minimum in R) for finite 
noise [ill Il7j . So there is an initial range of a- values 
where by increasing the noise the spike emission is facil- 
itated (T a reduces accordingly to Eq. ©) and this ren- 
ders the response of the neuron more and more regular. 
On the other side at high er-values the response becomes 
again more irregular, because the noise induces firing 
even during the relative refractory period thus modifying 
even the duration T e of the spike itself. To summarize, 
the activation process responsible for firing is gradually 
substituted by a diffusive process with drift and at the 
transition from one mechanism to the other a regulariza- 
tion of the output signal is observed. As shown in Fig.|I| 
for different values of I, both the coefficient of variation 
R and the signal correlation time r c were able to identify 
clearly the first kind of coherence in the silent regime. 
We can also see that t c increases and R decreases when 
I — * Isn, i-e., the coherence of the emitted spike trains 
increases by approaching the repetetive firing regime. 
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FIG. 4: Coherence of the emitted spike trains for increasing 
values of I in the silent regime, (a) r c (from bottom to top 
7=-l,0,l,2,3,4,5 /ttA/cm 2 ). (b) R (from top to bottom 7=- 
1,0,1,2,3,4,5 /iA/cm 2 ). 

The second kind of resonance could only be detected by 
t c . In fact this property is not related to spike emission 
(suprathreshold dynamics) but to the behaviour of the 
signal below the firing threshold. In Fig[S^ the behaviour 
of r c for I — AfiA/cm 2 is shown in a wider range of 
noise with respect to Fig. 0Jl In this range r c exhibits 
indeed two maxima: the first and higher maximum at low 
noise (a « 3) is related to the coherence of subthreshold 
oscillations, while the second one at intermediate noise 
(a s» 33) is due to the CR just discussed above. 

Let us now gain deeper insight about the origin of these 
two maxima. For extremely low noise (er < 3) the neuron 
rarely fires, i.e., T a — > oo; increasing the noise the sub- 
threshold oscillations generated by the input kicks are 
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FIG. 5: / = 4fj,A/cm 2 : (a) correlation time t c (empty circles) 
as a function of the standard deviation of the input noise a: 
(b) autocorrelation function C(t) for three increasing noise 
values a = 2.9 (continuous line), 4.5 (dotted line) and 9.7 
(dashed line); the vertical dot-dashed lines indicates the pe- 
riod Tj and the period Tf p for a — 9.7. 



more and more correlated. For a > 3 the occurrence 
of rare spikes tends to decorrelate the signal leading to 
a decrease of r c . With a further increase of the noise 
V{t) becomes essentially a sequence of spikes and C(t) 
represents the correlation between them; a second peak 
appears in r c indicating the maximal coherence of the 
spike sequence. 

To investigate the transition from one to the other 
behaviour let us analyze the autocorrelation function of 
V(t) for increasing value of a, as shown in Fig. [SJd. For 
a w 3 the autocorrelation function C(t) oscillates with 
a period w Ti (see Sec RD) . while at a » 9.7 the max- 
ima of C(t) correspond to multiples of Tf p , i.e., the first 
peak of Pisi{t). In between these two values there is a 
transition where the course of V(t), initially consisting 
of mere subthreshold oscillations, is dominated by the 
spikes. The transition is observed for a w 4.5, where in 
correspondence to the first oscillation C(t) reveals two 
maxima, one located at t « T; and another at t w Tf p . 



D. Conclusions 

We have presented an analysis of the response of the 
silent HH model in the high-input regime. In this regime 
the HH neuron displays a large variety of dynamical be- 
haviors and its response is completely determined by the 
average and the variance of the stochastic input. Our 
main result is the coexistence of two different coherence 
resonances: the first one, at intermediate noise levels, as- 
sociated with the regularization of the spikes emitted by 
the neuron; the second one, at very low noise, related 
to the coherence of subthreshold oscillations around the 
fixed point. The first one can be revealed using both the 
ISI coefficent of variation R and the autocorrelation time 
of the signal r c . The second type of resonance, observed 
for the first time, is not related to spike emission and can 
thus only be detected by means of r c . 

Moreover, we have examined the behaviour of the HH 
neuron for low and high input noise variance. The rich- 
ness of the HH dynamics is particularly pronounced for 
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low input noise where we have found the coexistence of 
a multimodal structure and an exponential tail in the 
Pisi(t). The peaks in the Pisi(t) suggests that the sys- 
tem, under the influence of stochastic inputs, can res- 
onate when forced with specific frequencies: the first peak 
is associated to frequencies in the 7-range [24| (namely, 
from 40 to 66 Hz for I £ [0 : 9] fj,A/cm 2 ), and the second 
one to lower frequencies (namely, from 30 to 37 Hz for 
the same interval of currents) 

In the silent regime the responsiveness of the single 
neuron is enhanced by stochastic stimulations, in partic- 
ular the regularity of the emitted spike trains is maximal 
in correspondence to an optimal noise level. Moreover, 
we expect that a population of such neuronal elements 
will have the capability to exhibit coherent and corre- 
lated activity over different time scales (mainly in the 



7 and /3-ranges 0), a property that is believed to be 
important for information encoding for cortical process- 
ing [23. Indeed it has been found 27] that for suffi- 
ciently strong synaptic coupling a globally coupled HH 
network subjected to stochastic inputs reveals a collec- 
tive synchronized rhythmic firing in a range of 40 — 60 
Hz, induced via CR. 

To proceed in the direction of more realistic situations 
the present analysis performed for a current-driven model 
should be extended by considering conductance-based 
synaptic inputs [2(1 Hfj. It is surely worth to address 
this issue in forthcoming studies, because the response 
of these two classes of models to noise fluctuations can 
sometimes be even opposite, as shown in |25j for the out- 
put firing rate. 
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